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1 INTRODUCTION 



ABSTRACT 

In the leptonic scenario for TeV emission from a few well-observed shell-type TeV 
supernova remnants (STTSNRs), very weak magnetic fields are inferred. If fast-mode 
waves are produced efficiently in the shock downstream, we show that they are viable 
agents for acceleration of relativistic electrons inferred from the observed spectra even 
in the subsonic phase, in spite that these waves are subject to strong damping by ther- 
mal background ions at small dissipation scales. Strong coUisionless non-relativistic 
astrophysical shocks are studied with the assumption of a constant Aflven speed in 
the downstream. The turbulence evolution is modeled with both the Kolmogorov and 
Kraichnan phenomenology. Processes determining the high-energy cutoff of nonther- 
mal electron distributions are examined. The Kraichnan models lead to a shallower 
high-energy cutoff of the electron distribution and require a lower downstream density 
than the Kolmogorov models to fit a given emission spectrum. With reasonable param- 
eters, the model explains observations of STTSNRs, including recent data obtained 
with the Fermi 7-ray telescope. More detailed studies of the turbulence generation and 
dissipation processes, supernova explosions and progenitors are warranted for better 
understanding the nature of supernova shocks. 
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The acceleration of cosmic rays up to ~ 10^^ eV has been attributed to super n ova explosions and TeV emission is ex- 
pected from the remnants i Ginzburg fc Ptuskin 19761 : Lagage fc Cesarskv 19831 : Reynolds 20081 : Butt et al. 20091 ) . The 
standard diffusive shock particle acceleration (PSA) model has been successful in explaining emissions from most supernova 
remnants jEichler Ill979cl : iBlandford fc Eichler Ill987l : Ixirk fc Duffv1ll999l : Izirakashvili fc Aharonian I boO?! : IVannoni et al 



2OO9I : iRevnolds Iboosh . Investigations of acceleration by a spec trum of turbulent p lasma waves, the so-called stochasti c par 



tide acceleration (SA), also have a long and resilient history JScott_&Chevalier_ 1975| 



Eilek 



197S 



Bvkov fc ToDtvein Ill983l: IC owsik fc Sarkar Ill984l: 



Petrosian fc Liu 



2004 : 



Cho fc Lazarian . 2006: Liu et al. 



Ptuskin 



Lacombe 



1977 



Achterberg 



19881 : iBvkov fc Toptygin Ill993l : lAtovan et al. 



197S 



200c 



2008a|). Most authors prefer the use of relativistic leptons to ac- 



count for the nonthermal r adio. X-rays, and TeV emissions from the remnants. The TeV emission has also been attributed t o 
energetic protons and ions l Aharonian et al. 20061 : Morlino et al. 20091 : Fang et al. 20091 : Zirakashvili fc Aharonian 20ld ). 

Although the DSA can naturally give a universal power-law energetic particle distribution with the spectral index de- 
termined by the shock compression ratio for a lineax model, it requires well-defined shock structure and efficient scatter of 
high-energy particles by small-scale turbulence ( Bell 19781 '). The initial acceleration of low-energy particles to a high e nough 
energy for the shock to be effective, the so-called injection problem, is likely due to the SA by turbulence. iLacombe I (|l977l 'l 
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showed that the S A by smaU-scale Alfven waves can be more efficient than the first order Fermi acceleration by shocks. 
Achterberg (|l979l ') derived approximate diffusion coefHcients and showed that stochastic interactions of particles with a 



plored for broader astrophysical applications (e.K-, Eichler 


1979a b; Cowsik & Sarkar 19841: Ball. Melrose. & Norman 1992: 


Miller. LaRosa & Moore 


1996: Schlickeiser & Miller 1998 


:|Petrosian & Liu 1 2004: Yan & Lazarian 1 20041: Fisk & Gloeckler 


2OO7I: Petrosian & Bvkov 


20081). The essential challenges 


to the SA model are a self-consistent treatment of the nonlinear 



turbulence spectral evolution and the requirement of t he same energy dependence of the acceleration and escape timescales 
for th e production of a power-law particle distribution ( Ball. Melrose, fc Norman 19921 : Becker. Le fc Dermer lliood : Ptuskin 



1988h . Rec ently Gibbsian theory has bee n generalized to account for power-law distributions in marginally stable Gibbsian 



equilibria ( Treumann fc Jaroschek 20081 '). It remains to be shown how the physical processes of the SA are related to the 
ordering pa rameter k of this statistics. Most previous studies assume isotropic magnetohydro dynamic (MHD) waves for the 
turbulence jLacombe Ill977l : Uchterberg II1979I : lEichler Ill979bl : iMiller. LaRosa fc Moore Ill996l'). The anisot r opy of the turbu- 
lence caused by cascade and damping proces s es has been consid ered recently I Goldreich fc Sridhar 1995 : Chandran 20031 : 



Yan fc Lazarian I I2OO4I : ICho fc Lazarian Ibood : Ijiang et al.1l2009l ') . In particular, iLiu et al. I l|2008al ) show that the scatter and 
acceleration rates of charged relativistic particles by fast-mode waves in a high-/3 plasma may be much higher than those 
given by the standard quasi-linear theory with an isotropic wave power spectrum. 

Over the past few years, detailed radio. X-ray, 7-ray, and TeV observations of a few shell-type TeV supernova remnants 
(STTSNRs) pose several challenges to the classical DSA mode l in the hadronic scenario, where the TeV emission is produce d 
through neutral pion decays induced by proton-proton scatter ijAharonian et al. I2OO6I . l2007l : iTanaka et ani2008l : iFunk II2OO9I ) . 
Besides requiring efficient amplification of the magnetic field in the upstream plasma and a good correlation between the 
magnetic field and background plasma density ( Plaga 20081 : Fang et al. 20091 : Morlino et al. 20091 ). the model also implies 
a cosmic ray energy of ~ 10^^ ergs for e ach remnant and a n electron acceleration efficiency more than 4 orders of magnitude 
lower than the ion acceleration efficiency (|Butt et al. 1120081 '). The high density of the upstream plasma in the model also implies 



signif icant thermal X-ray emission from the downstream, which may exceed the observed upper limit (jCassam-Chenai' et al. 



.2004h. A very hard proton spectrum is also required to fit the 7-ray spectrum obtained recently with the Fermi 7-ray telescope 
jFunklboogf i. Although these kinds of remnants ma y be atypical, detailed modeling can still have profound implications for 
our understanding of supernova shocks ( Butt et al. 20091 : Zirakashvili fc Aharonian 2O10l ). 

In a previous paper, we showed that the SA of el ectrons by turbulent plasma waves in the shock downstream might 

2008al ). Turbulence is expected giv en that the size of these remnants ar e 



naturally explain these observations (jLiu et al 



many orders of magnitude larger than the dissipation scale of the ion inertial length ( Dickel et al. 1991 : Jiang et al. 2009l l 



The DSA rnodel proposes that particle acceleration occurs directly and predominantly at the ion inertial length or gyro-radius 
jRevnolds 2008li. This req uires the absence of instabilities over a large range of spatial scales, which is highly idealized. 



Giacalone fc Jokipii I (|2007l ) showed that density fluctuations in the upstream can be amplified significantly by shock waves. 



resulting in strong turbulence in the downstream. Magnetized turbulence appears to be a more generic and natural energy 
dissipation channel than the short-length-scale shock fronts (SF). 

With the leptonic model, the TeV e mission is mostly pr oduced by the inverse Comptonization of the cosmic microwave 



background radiation by TeV electrons I Porter et al. 20061 ) . The magnetic field required to reproduce the observed X-ray 



flux by the same TeV electrons through the synchrotron process implies a spectral cutoff in the hard X-ray band, which is 
in agreement with observations. The SA model also requires a much lower gas density than the DSA model, which not only 
explains the lack of thermal X-ray emission from the shell of the remnants, but also reduces the energetics of th e supernovae 



The required turbulerice ge neration scale is comparable to the size of the observed X-ray fllaments as well ijDickel et al. 



energy electrons (|Liu et al 



I991I : lUchivama et al. 2007 ) . The fast variability of small X-ray features may be attributed to rapid spatial diffusion of high 



2008a|). As we will show in this paper, the agreement between preliminary results from Fermi 



observations and the model prediction is also impressive |Funk II2OO9! ) 



\ ' 1 1 '■ — I ' ' 

The SA by fast-mode waves has b e en studied by several a,uthor s (IBvkov fc Toptvg in 1983: IPtuskin 1 1988: Bvk oy fc Toptygin 



gu 

1993 : iMiller. LaRosa fc Moore] Il996l : ISchlickeiser fc Miller |[l998l : lYan fc Lazarian 1.2004 : ,Cho fc Lazarian 20061 : IlIu et al. \ 
20061 ). Both resonant and nonresonant interacti ons have been considered. I n this pap e r, we consider the nonresonant accel- 



eration by compressional waves flrst studied bv iBvkov fc Toptvgin I (| 19831 '): iPtuskin I (|l988l ). Compared with these original 
studies, our model has several distinct features: 1) most of the dissipated fast-mode turbulence energy is absorbed by thermal 
background ions; 2) the residual fast-mode waves in the dissipation scales propagate along local magnetic fields and preferen- 
tially accelerate electrons in the background plasma; 3) the plasma physics processes in the dissipation range in principle may 
lead to a self-consistent treatment of the electron injection process at low energies; 4) the scatter mean-free-path of relativistic 
electrons, which is a free parameter in most of the previous studies, is determined by the characteristic length of the magnetic 
field, which in a high-/3 plasma is reduced by strong turbulence motions significantly; 5) the high-energy cutoff of the particle 
distribution is tied to the characteristic length of the magnetic field. 

In this paper, we first discuss the SA by decaying turbulence in general and show that fast-mode waves may account for 
observations of a few STTSNRs (Section [2]). In Section [S] we present the structure of the downstream turbulence with both 
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the Kolmogorov and Kraichnan phenomenology for the turbulence cascade. The transit-time damping (TTD) by the thermal 
background particles of compressional fast-mode turbulence is considered in the dissipation range. The stochastic acceleration 
of relativistic electrons by fast-mode wave turbulence in the subsonic phase is discussed in Section [4l where physical processes 
determine the acceleration of the highest energy electrons are discussed. The models are applied to the well-observed TeV 
SNR RX J1713. 7-3946 in Section [5l Although the inferred plasma density is much lower than that in the hadronic scenario, 



some models still have significantly higher densities than that derived from XMM-Newton observations (jCassam-Chenai' et al. 



200J). In Section [6l we discuss how the density may be further reduced by considering the turbulence generation processes 



and first-order Fermi acceleration in the supersonic phase. Conclusions are drawn in Section [T] 



2 GENERAL CONSTRAINTS ON THE PARTICLE ACCELERATION 

According to the classical Fermi mechanism ( Fermi Ill949l ), the particle acceleration rate is determined primarily by the scatter 
mean- free-path / and the velocity of the scatter agents u ( Blandford fc Eichler 1987 : Ball. Melrose, fc Norman 19921 ). Very 
general const raints can be ob tained on the nature of these processes when applying this mechanism to specific observations. 
For example, lEichler I (|l979al ) showed that the acceleration process must be selective in the sense that only a fraction of the 
background particles are accelerated to very high-energies. Otherwise, these stochastic interactions likely lead to plasma heating 
instead of a very broad energy distribution of accelerated partic les as frequently observed in dynamically evolving coUisionless 
astrophysical plasmas. For particle acceleration in solar flares, Eichler (|l979bl ) argued that this selective acceleration could 
be achieved in the energy domain (i.e., the frequency domain for the waves) through cyclotron resonances of particles with 
a spectrum of cascading plasma waves. In the PSA model, the pa rticle acceleration at low energies and the shock structure 
determine the efficiency of different particle species (jEichler 1 ll979al ). In the presence of a magnetic field, selective energization 
may also be realized in the do main of the particle pitch-angle and/or wave direction angle with respect to the magnetic field 



Beresnvak fc Lazarian 2008h 



Further constraints can be put on the SA by a spectrum of turbulence. Given the small gyro-radii of charged particles in 
magnetized astrophysical plasmas, charged particles couple strongly through the magnetic field. As a result, the turbulence 
responsible for the SA will decay as the energy carrying plasma being carried away from the source region of the turbulence 
by large scale flows and /or magnetic fields for a high and/or low f5 plasma, respectively. This is the case for the SA in a shock 
downstream with a high value of the plasma /3, where the turbulence is generated at the SF and its intensity decreases as the 



flow moves away from the SF (jLiu et al 



In the Kolmogorov phenomenology for the turbulence cascade (jKolmoeorov 

given by 

Q = Cipu^/L , 



20083). We next discuss co nstraints on such a particle acceleration scenario. 

" l94ll ). the free energy dissipation rate is 



(1) 



where Ci is a dimensionless constant, p is the mass density, and u and L are the eddy speed and the turbulence generation 
scale, respectively. The eddy turnover speed and time at smaller scales are given respectively by 



47rW(fc)fc^ oc k~^^^ , 



Vedd{k) 
Tedd{k) 

where 

W{k) = (?iV47r)i-'/'fc-"''' = (4-Ky\Q/Cipf^^k-^^^^ oc fc-"/' 



(2) 
(3) 

(4) 



is the isotropic turbulence power spectrum, fc = l/l is the wave number and I is the eddy size. From the three-dimension 
Kolmogorov constant C ~ 1.62 ( Yeung fc Zhou 1997), we obtain Ci = C"'/^ = 0.485. At the turbulence generation 
scale L = 1/fcm, v^dd = u, Q = C\p\(A-KW)^k^^]^^'^ = C\pv1^^{k) /r^ddik) , and the total turbulence energy is given by 
j VK(fc)47rfc^dfc — (3/2)u^ . The turbulence decay time is therefore given by Td = df/dln(u) = ireddikm) /Ci , where t 
indicates the time, i.e., eddies decay after making 3/(27rCi) ~ 1 turn. 

We are interested in the accelera tion of particles throu gh scattering randomly with heavy sc atter centers with the corre- 
sponding acceleration time given by jBlandford fc Eichler ||1987. : . Ball. Melrose, fc Norman Ill993 ') 



= [3W /Wedd(fc)] : 



where 



= (kv) ^ = l/v 



(5) 



(6) 



is the scatter time, v is the particle speed, and we have assumed that the scatter mean free-path is equal to I. For the above 
isotropic Kolmogorov turbulence spectrum, Tac(k) — 3v / (AtvW k"^) oc fc"^/' . To have significant stochastic particle acceleration, 
the acceleration time Tac{k) should be shorter than the turbulence decay time Td, which implies v? > Civvedd{k) . So, in general, 
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the SA is more efficient at smaller scales. The onset scale of the SA is given by kc = {C\v/u)^km , where Tac = t^. Therefore, 
to produce energetic particles with a speed of u by a Kolmogorov spectrum of scatter centers, the turbulence must have a 
dynamical range greater than 

Dkoi = {Civ/uf . (7) 
This acceleration process corresponds to the acceler ation by incompre s sional motions studied by Bvkov fc Toptvein ( 19831 ). 



In the Iroshnikov — Kraichnan phenomenology (jlroshnikov 



19631 : 



Kraichnan 



19651 : 



Jiang et al. 



20oi), the turbulence 



cascade rate is suppressed by the wave propagation effect by a factor of vp/vedd, where vf (independ of k) is the wave group 
speed: 



CipV^ad/ivFTedd) = Q ■ 

Then we have 



Q 


= C-ipu' /(Lvf) , 


W{k) 


= {u^/A-K)k]i^k-''^^ 


Vedd 




Tedd 






= (3«/^2)(fc„fc)-^/2 



and the turbulence decay time is given by 

Td = 'iT^M{km)VF / {Clu) , 



(8) 

(9) 
(10) 

(11) 
(12) 
(13) 

(14) 



where the wave speed vf 2> ii. To have significant acceleration through scatter with the eddies, the dynamical range of the 
turbulence must be greater than 



DiK = {Civ/vpf 



(15) 



which is much less than DkoI = {Civ/u)^. 

The resonant interactions of particles with waves may be more effective in accelerating particles in this case than the 
interactions with eddies. For a wave phase speed vf independent of k and a wave spectrum of W{k) = ii|n(47rfcm)~^(fe/fcm)~*, 
where S is the wave spectral index, the standard quasi-linear theory gives a scatter time of Tsc — {L /v)(rgkmY~^ — 
\v%/vedd{r'^^)'^]rg/v for particles with a gyro-radius of rg = k^^ . And the acceleration time is given by Tac — 3tscV^/vf — 
{3v'^ /vedd{ko)'^)/{vko), which is essentially the same as equations ((5]) and ((6]). Although the acceleration rate depends on 
the phase speed of these waves, the scatter rate is proportional to the intensity of waves in resonance with the particles. So 
the acceleration is not enhanced by these kinds of r esonant interacti o ns. In the presence of m agnetic fields, the stochastic 
acceleration of charged particles can be very efficient (jchandran Il2003l '). lcho fc Lazarian I (|2006l ) reviewed several mechanisms 
by MHD waves, which can accelerate particles within the turnover time of large scale eddies for appropriate particle spatial 
diffusion coefficients. Some of these mechanisms likely result in plasma heating instead of particle acceleration as the particle 
energization processes are not selective. The corresponding waves are also likely subject to efficient damping by the thermal 
background particles. There are other mechanisms, which preferential ly accelerate en ergetic particles. 

For isotropic acoustic wave turbulence with a wave phase speed vf, Ptuskin 1 1988h showed that the st ochastic acceleration 
timescale of energetic particles with a spatial diffusion coefficient D = TscV'^ /3 > VF/kd is given by (see also lBvkov fc Toptvein 



[3 + aS){Dk^/vFf-']D/u^ 



(16) 



where 3u^/2 is the overall wave intensity, km indicates the wave generation scale, and below the turbulence dissipation scale 
Id = 1/kd the wave intensity is negligible. ^(5) depends on S and is on the order of unity (e.g ., ^(11/3) 3.0). The results for 
fast a nd slow diffusion have been combined here approximately to give a unified expression ( Ptuskin 19881 : Cho fc Lazarian 
2006[ ). For 5 < 4, Tac increases monotonically with D. To have significant acceleration, the minimum acceleration timescale 
= TscV^ /y? must be less than Td — ^Lvf / {C\v?). Then the scatter mean-free-path of the particles TscV must be shorter 



Tacuiii 

than 



1/2 



L/Dai = iLvF/{C^v) = 3L/Dy^ 



(17) 



It is also possible that the dynamics of the turbulence cascade is not affected by the wave propagation, which only enhances 
the particle acceleration rate. Then Td = 3I//(d?i), the dynamical range required for the SA to be significant is given by 

Dak = C^v/iiu) = Dj/'i/S . (18) 



Several STTSNRs have been observed extensively in the radio. X-ray, and TeV bands. X-ray observations with the 
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Chandra, XMM-Newton, and Suzaku, and TeV ob servations with the HESS have made several surprising discoverie s that 
chaUenge the DSA model in the hadronic scenario ( Aharonian et al. 20061 . 2007 : Liu et al. 2008a : Tanaka et al. 20081 ). The 
leptonic scenario, on the other hand, is relative simple except that the ele ctron acceleration mechanism needs to be addressed 



I Zirakashvili fc Aharonian 



2007: Liu et al. 2008a: Vannoni et al 



20091 ). 



The TEV SNR RX J1713.7-3946 is about Tufe 



1600 years old (|Wang et al. I ll997l l with a radius of 7? ~ 10 pc and a distance of D ~ 1 kpc. By fitting its broadband spectrum 
with an electron distribution of / oc 7"'' exp — (7/7c)^'''^ (the dashed lines in Fig. [T]!, we find that p — 1.85, B — 12.0 /iG, 
'ycTTieC? = 3.68 TeV, and the total energy of relativistic electrons with the Lorentz factor 7 > 1800 = 3.92 x 10'*^ erg. 

The X-ray emitting electrons have a gyro-radius of ~ 10^^ cm, which should be shorter than the particle scatter mean- 
free-path. To produce these electrons through the SA, the turbulence must be generated on scales greater than DKoiTg, DjKVg, 
DarTq, and DAiVg for the non-resonant Kolmogorov, Kraichnan, acoustic Kolmogorov, and Kraichnan phenomenology, 
respectively. For STTSNRs, u ^ vp ^ 0.01c, DkoiTq ~ 10 kpc, which is much greater than the radii of the remnants. 
The SA by eddies with a Kolmogorov spectrum is therefore insignificant. The standard quasi-linear theory also predicts 
negligible SA. DjKTg ^ 30 pc, which is also too thick. DAiTg ~ DarTq ~ 0.03 pc, which is much greater than the ion 



inert i al length and may result fr o m Kelvin-Helmholtz instabilities or cosmic ray drifting upstream ( Bell 19781 : Micono et al 



I999I : ICiacalone fc Jokipii |[2007l : iNiemiec et ahlboO^ ). Therefore if relativistic electrons from the STTSNRs ar e accelerated 



through the SA, they may be energized by high-speed compressional plasma waves. Previously lEichler I (|l979al ) showed that 
the non-selective acceleration as given by equations ([5} and (|13p leads to plasma heating instead of particle acceleration. 
Here we give a different argument against these acceleration processes. The acceleration timescale in the following is given by 
equation (I16|l with D depending on the wave spectrum and particle energy. 

The DSA models are usually favored over the SA models for two main reasons: 1) the DSA can naturally produce a 
power-law particle distribution; 2) the DSA corresponds to a first-order Fermi mechanism and is presumably more efficient 
than the SA, which corresponds to a second-order Fermi mechanism. While how shocks produce power-law high-energy particle 
distributions is relatively well understood, the second reason appears to be a misconception. It is true that the ratios of the 
acceleration and scatter rates are proportional to the first and second power of the speed ratio of the scatter agent and the 
particle for the DSA and SA, respectively. But in the DSA models, there are two scatter processes: scatter of particles by 
turbulence, which causes the particle diffusion, and the particle crossing of the SF due to this diffusion. The SF crossing rate 
is ab out a factor oi v/u lower than the particle scatter rate by the turbulence for a shock speed of ~ w i Lagage fc Cesarskv 
19831 ). As a result, the acceleration rates of the DSA and SA models are both proportional to /Iv. Therefore the SA is 



not necessarily less efficient than the DSA. The SA at small scales can also be enhanced by high speed kinetic plasma waves 



i Prvadko fc Petrosian Ill997l : IPetrosian fc Liu II2OO4I ) 



Of course, the speed of the scatter agent accessible to a particle may 
depend on the scatter mean-free-path I in the SA. The stochastic acceleration rate can be very low if decreases quickly with 
the decrease of I. For the DSA and SA by large-scale acoustic waves in the shock downstream, the speed of the scatter agent is 
proportional to the shock speed. One just needs to reduce the scatter mean-free-path / to enhance the acceleration rate. As we 
will show in Sections [3] and [H considering the competition of isotropic cascade with anisotropic damping of fast-mode waves 
through the transit-time-damping process, the waves in the dissipation range propagate along local magnetic field lines with 
a spectrum index of 2. In a high-/3 plasma, the damping is mostly caused by ions. The residual parallel propagating waves 
however preferentially accelerate electrons. Enegetic particles are also subject to efficient acceleration by large-scale fast-mode 
(acoustic) waves with energy independent acceleration and scatter rates, giving rise to power-law particle distributions in the 
steady-state. At even higher energies, the gyro-radius of the particles exceeds the characteristic length of the magnetic field, 
which also corresponds to the dissipation scale, the spatial diffusion coefficient D increases quickly with energy. However, 
due to interactions with large-scale turbulence, both the acceleration and spatial diffusion timescales vary gradually with D, 
implying a gradual high-energy cutoff. 



3 SHOCK STRUCTURE AND DAMPING OF FAST-MODE WAVES IN THE DOWNSTREAM 

We consider the relatively simple case, where both the thermal pressure and magnetic field are negligible in the upstream 
and the shock normal is parallel to the plasma fiow. The mass, momentum, and energy fiuxes are given, respectively, by 
pVo, P + pVo , and Vo{E + P + pVo /2), where Vo, P, and E are the speed, pressure, and energy density of the plasma flow, 
respectively. Based on the parameters inferred from the leptonic scenario for the TeV emission from STTSNRs shown in 
Section [2] the plasma /3 in the shoc k downstream is l ikely high and the fast-mode wave damping by thermal background is 



then dominated by protons and ions I Liu et al. 2008a ). The Alfven speed given by va ~ {B^ /A-Kp)^^^ therefore must be much 
smaller than the turbulence speed u near the shock front, where B is the magnetic field intensity. For strong non-relativistic 
shocks with the shock frame upstream speed U much higher than the speed of the parallel propagating fast mode waves in 
the upstream vp = {va + 5u|/3)^''^, where = Pg/p is the isothermal sound speed and Pg is the thermal pressure of the gas, 
Vo — U and pU^ ^ P ^ E in the upstream. In the downstream, the pressure and energy density have contributions from the 
thermal gas and turbulence: P = p{v's + ii^) and E — p{3vg/2 + 3u^/2), where we have assumed that the turbulence behaves 
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0123456789 10 0123456789 10 

X/L X/L 



Figure 1. Evolution of tlie eddy speed u and the speed of parallel propagating fast mode waves vp in the downstream. The thick and 
thin lines are for the Kraichnan with = 0.036C/ (left) and Kolmogorov with = O.OlSf/ (right) phenomenology, respectively. 



as an ideal gas and ignored possible dynamical effects of the wave propagation. Then we have (|Tidman fc Kraillll97ll ) 



puU = pdU/A, 
PuU^/2 = pd(7((7Vl6 + 54 +5u'')/8, 



(19) 
(20) 
(21) 



where the subscripts u and d denote the upstream and downstream, respectively, and we have ignored the effects of the 
electromagnetic fields and the thermal energy and pressure in the upstream. Then we have 



(22) 



This is slightly different from that given by iLiu et al. I (j2008al l. where we assumed that the pressure and enthalpy of the 
turbulent magnetic field are given by pdVA/'i and 5pdU^/6, respectively. 

The shock structure can be complicated due to the present of turbulence. We assume that the turbulence is generated 
isotropically and has a generation scale of L, which does not change in the downstream. The speeds vs, va, and u therefore 
should be considered as averaged qua ntities on the scale L. va depends on the upstream conditions and /or the dynamo 
process of magnetic field amplification ( Lucek fc Bell 2OO0I : Cho fc Vishniac 200Cll : Niemiec et al. 2008h . Here we assume it 
a constant in the downstream. One can t hen quantify the evolution of other speeds in the downstream. 

For the Kolmogorov phenomenology ( Zhou fc MatthaeusI 1990l ) , 



iApv? 
2At 



■jUAujx) 
8dx 



Ciu{xf 



(23) 



Near the SF, we denote the isothermal sound speed and Afiven speed by vso and vao, respectively. Then the eddy speed at 
the SF is given by a^^'^U/A with a = 3 — 16«|o/t^^- Integrate equation (|23|) from the SF (x — 0) to downstream {x > 0), we 
then have 



u{x) 




1 


U 4Cia;/3L + 4/ai/2 ' 


vs{x) 


' 3 


1 


U 


16 ' 


16(Cix/3L + a-i/2)2 


vf{x) 


' 5 


5 


U 


16 ' 


48 (Cia:/3L + a-i/2)2 



1/2 



+ 



[/2 



1 1/2 



(24) 
(25) 
(26) 



As mentioned in Section (2] to produce the observed X-ray emitting electrons in the STTSNRs through the SA processes, 
fast-mode waves needs to be excited efficiently. The MHD wave period is given by rpik) = 2-e jvph. Then the tr ansition from 
the K olmogorov to Kraichnan phenomenology occurs at the scale, where Tpikt) ~ 2nTcdd{kt) or vp = Ucdd(fct) (jJiang et al 
2OO9I) . We then have 



kt = (u/vf) k„ 



(27) 
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For kt > k > km, the turbulence spectrum is Kolmogorov like given by equation Q. For k > kt > km, the turbulence 
spectrum in the inertial range is given by 



W{k) 



_1_ 

47r 



vj. u ' k 

in 



2 1.2/3,-11/3 



u k 



for IK , 
for Kol . 



(28) 



Although the turbulence energy exceeds (3/2)u^ when the wave propagation effect is considered, we still assume that the 
enthalpy of the turbulence is given by (5/2) u'^ for vf < u so that equation p2|l and the above solutions for the speed profiles 
remain valid. 

In the subsonic phase with vf > u, 

kli^k-'^!'^ for IK, 

for Kol . 



(29) 



and 

3Udu{xf 
8da; 



Cizt(a)-^ 
Lvf 



for IK , 



(30) 



where from equation (|22|) one has vf = [5(7^/16 + t"^ — 5u^ (x) /3] . Equation (O can be solved numerically to get the 



speed profiles in the subsonic phase. Figure [T] shows the vf and u profiles with va = vao 
thick and thin lines are for the Kraichnan and Kolmogorov phenomenology, respectively. 

In summary, for k > max (km, kt), the turbulence spectrum in the inertial range is given by 



vso U in the downstream. The 



W{k) 



1 

An 



3/2 • /■ 1/2 

u ' mm[vp , u 



l/2^I.l/2j, -7/2 



for IK , 
for Kol , 



(31) 



and for kt > k > km, W{k) = u'^ (Atv)-'^ kK^ k-'^^^'^ . 

The transit-time damping (TTD) of compressional fast-mode waves starts at the characteristic length of the magnetic 
field Id = l/kd , where the Alfven speed is comparable to the eddy speed, i.e.. 



Vl = i7TW{kd)k^d 



dn{v]/\u^^^)u^^\U'k-'/^ for IK, 

\km/kif'^ for Kol. 



(32) 



At even larger scales, the vortex motions produce random magnetic fields comparable with the mean field reducing the scatter 
mean-free-path of charged background particles to Id or even shorter scales. This trapping of charged particles within a scale 
of Id prevents the TTD on scales above Id- The inco mpressional modes are not subject to the TTD and have different spectra 
( Goldreich fc Sridhar |[l995l : Icho fc Lazarian Ibood ') . ' 'n what follows, we only consider the compressional (fast magnetosonic) 
wave modes. Then we have 



min(iiF, u)/vA 
1 



for IK , 
for Kol . 



(33) 



Since vf > va, one has kt < kd, and damping is negligible in the regime, where kt > k > km- 

For a fully ionised hydrogen plasma with isotro pic particle distributions, which is reasonable in the absence of strong 
large-scale magnetic fields, the TTD rate is given by jStix II 19621 : iPetrosian et al. 1120061) 



AT(6i,fc) = 



(27rfcB)^^^fcsin^6' 
2{m,e -f rUp) cos 9 



1/2, 



II -I- (Tpmp) ' e 



(34) 



where fee, Te, Tp, me, rup, 6, ui, and k\\ — kcosO are the Boltzmann constant, electron and proton temperatures and 
masses, angle between the wave vector and mean magnetic field, wave frequency, and parallel component of the wave vector, 
respectively. The first and second terms in the brackets on the right hand side correspond to damping by electrons and protons, 
respectively. For weakly magnetized plasmas with va < vs, proton damping always dominates the TTD for LtP' jk'^ 



kBTp/mp-ii If VA does not change dramatically in the downstream, the continuous heating of background particles through 
the TTD processes makes Tp {mp/me)Te since the heating rates are proportional to (mT)^^'^, where m and T represent the 
mass and temperature of the particles, respectively. We see that parallel propagating waves (with sm6 — 0) are not subject 
to the TTD processes and can accelerate some particles to relativistic energies through resonant interactions. Obliquely 
propagating waves axe damped efficiently by the background particles. Although the damping rates for waves propagating 



1 For Te = Tp, the electron damping term dominates when va > 1.9 vg. 
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Figure 2. The angular-integrated compressional mode turbulence spectra v'^^^ik) at several locations in the downstream indicated in 
the legend for the speed profiles in Figure [T] In the dissipation range, it corresponds to fcVV instead of f^^^- The onset of the TTD at 
causes the discontinuity for the former. The left panel is for the Kraichnan phenomenology. The Kolmogorov, Kraichnan, and damping 
ranges are indicated for the supersonic phase spectrum with x = O.OIOL. At the other locations, the turbulence is subsonic and there are 
only Kraichnan and damping ranges. The right panel is for the Kolmogorov phenomenology. 



nearly perp endicular to the maRn etic field (cos S ~ 0) are also low, these waves are subject to damping by magnetic field 
wandering ( Petrosian et al. 2006^). 

The turbulence po wer spectrum cuts ofT sharply when the damping rate becomes comparable to the turbulence cascade 
rate ( Jiang et al. 20091 ) 

^7dd/(TF ' + -Tedd) - '^IdW for IK , , . 

1 forKol. ^^^> 

One can define a critical propagation angle 9c{k), where AT{Oc,k) = r(fc). Equations (|3ip and (|34p then give 



( 1 ^ / vlkl/'/i2-^/'.^f\sv,k^/') for IK , 

■''''PI 24cos2eJ - \ „^fcyV(2-'/V/2„sfci/3) forKol. ^ ' 



where the electron damping term has been ignored. The angular-integrated turbulence spectrum in the dissipation range is 
therefore given by 

el{k) { u^/^min{vy\u^'^)k]i\-^'^ for IK, 
W{k) ^ ^1 forKol. 

exp(5/6)«ifcd,_2/ va/vf for IK , 



2i/27ri/2„g 1^ 1 forKol. 

Interestingly the turbulence spectrum is inversely proportional to k^ in both scenarios. The angular-integrated turbulence 
spectra J ty(fe)27rfc^dcos6' for the velocity profiles in Figure [T] at several locations in the downstream are shown in Figure (2] 
The discontinuities of the angular-integrated turbulence spectra are caused by the abrupt onset of thermal damping at the 
characteristic length Id of the magnetic field. Obliquely propagating fast-mode waves are absorbed by the thermal background 
ions at this scale. 



4 STOCHASTIC ELECTRON ACCELERATION BY FAST-MODE WAVES IN THE DOWNSTREAM 



In a magnetized plasma, fast -mode waves are likely the agent responsible for efficient SA of electrons (|Bvkov fc Toptygin 
19831 . I1994. IChandran M2003|). The resonant interactions of particles with fast-mode waves have been studied by several 
authors jMiller. LaRosa fc Moore 19961 : Schlickeiser fc Miller 19981 : Petrosian fc Liu 2004). In these studies, the authors 
prescribed the wave spectrum with several parameters and calculated the corresponding Fokker-Pl anck coefficients. A self - 
consistent treatment of the turbulence spectral evolutio n in the dissipation range was presented by Yan fc Lazarian ( 200d ) 



for a variety of astrophysical plasmas. Liu et al. ( 2006h considered the damping of fast-mode waves in a low-/? plasma and 



the application of the SA of relativistic protons in magnetic field dominated funnels derived from general relativistic MHD 
simulations of non-radiative accretion fiow around black holes. Section [3] shows that only fast-mode waves propagating parallel 
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to the local magnetic field can survive the TTD by the thermal ions in the high-/? downstream plasma. These waves are right- 
handed polarized and resonate pr eferentially with electron s in the thermal background and therefore may selectively accelerate 
electrons to relativistic energies IjPetrosian fc Liu 1120041 ). We next explore the SA of electrons in the shock downstream in 
this scenario. 

We assume that the turbulence is isotropic from the turbulence generation scale L down to the dissipation scale Id ~ 1/kd, 
which depends on the turbulence speed u (Eq. [33]). Fast-mode waves are excited at the scale If = 1/ max (^kt J km ) — 
L/ max (1, and are also isotropic down to the dissipation scale. Due to the TTD, the energy density of fast-mode 
waves in the dissipation range is less than the magnetic field energy density. The resonant scatter rate of energetic particles 
by MIfD waves is therefore smaller than v/ld- Ifowever, for this high-/? plasma, particles with a gyro-radius r-g less than the 
characteristic length of the magnetic field Id draft along magnetic field lines with a scatter mean- free-path I = vt^c < Id- The 
scatter mean-free-path of particles with r-g ~ Id should be comparable to Vd — Id since particles with even higher energies 
scatter with the magnetic field randomly instead of performing gyro-motions. This efficient scatter of low-energy particles 
by the turbulence magnetic field results from bending of magnetic field lines by strong turbulence in the inertial range 
beyond Id- The resonant wave-particle interactions in the dissipation range give a much lower scatter rate due to the strong 
damping of obliquely propagating fast-mode waves by thermal ions. Although this scatter through the particle gyro-motion 
and chaotic magnetic field struture is not caused by resonances with fast-mode waves, it dete rmines the spa tial diffusion 
coefficient v'^Tsc/S — vld/S, which plays essential roles in the SA by large-sc ale fast-mode wave s dPtuskin 198S ). 

The SA in the supersonic phase with u > vf is not well understood I Achterberg 199G| : Bvkov fc Toptvgin 19931 ). If 
we assume the turbulence speed u as the characteristic speed of the scatter agents and a scatter rate of Id/c for relativistic 
particles with r-g < Id, where c is the speed of relativistic particles, the acceleration rate will be comparable to that of the 
DSA. We also note that the turbulence speed u is higher than vf in a narrow region near the SF (Fig. [TJ corresponding 
to the turbulence generation. Such an SA will be difficult to distinguish from the first-order DSA. Further downstream, the 
turbulence speed is lower than vf- We will ignore the particle acceleration in the supersonic phase a nd only consid er the 
acceleration by acoustic waves in the subsonic phase, where Equation (|16|) is approximately applicable ( Ptuskin 19881 ). 



We assume Tsc = Id/v = Id/c in the following for relativistic particles with r-g < Id- Particles with even higher energies 



have a scatter ti me Tsc > rg/c. The corresponding spatial diffusion coefficient is given by _D = TscC /3 ( Lagage fc Cesarskv 



19831 ). Following Ptuskin ] ( 19881 ) . we have the acceleration timescale of these particles by a spectrum of fast-mode waves given 
by 



~9~ 



dk- 



;2 +D2fc2 



(38) 



The fast-mode turbulence also enhances the spatial diffusion of these particles. The inc ompressional modes are more efficient 
in enhancing the spatial diffusion than compressional modes (IBvkov fc Toptvgin 1993 ) . To include these effects an d partially 
take into the effects of energy loss due to adiabatic expansion (|Bvkov fc Toptvgin Ill98i : ICowsik fc Sarkar Ill984l ). we adopt 
an effective diffusion coefficient: 

D, = D + xuL, (39) 

where X is a dimensionless parameter. The escape time of relativistic electrons from the acceleration region is then given by 

T,sc^{kl,D,)-^ =L^/D.. (40) 

Physically, x needs to be less than 1, which corresponds to maximum diffusion caused by turbulence. However, considering 
the possible presence of incompressional modes and energy loss due to adiabatic expansion, the acceleration timescale by 
compressional modes will increase since the incompressional modes will carry part of the turbulence energy. The effect of 
this increase of acceleration timescale on the electron distribution can be partially taken into account by reducing the escape 
timescale, i.e., by increasing the effective spatial diffusion coefficient D,. In what follows, we will treat x as a free parameter 
with X > 1 indicating the presence of incompressional modes and reduction of acceleration by fast-mode waves. The presence 
of incompressional modes will also bring the transonic point closer to the SF. The acceleration by compressible modes in the 
subsonic phase can sti ll be efficient. The details of these processes will depend on the coupling between compressible and 
incompressible modes ( Petrosian fc Bvkov 20081 ) and is beyond the scope of this paper. 

The dependence of acceleration and escape timescales on D at the transonic point in the downstream for typical conditions 
of SNR RX J1713. 7-3946 are shown in Figure[3l The approximate acceleration timescale given by Equation (|16|) and the escape 
timescale due to the diffusion coeflicient D alone, (fc^_D)~^, are indicated with the thin lines. At high values of D, the slightly 
high discrepancy in the exact and approximate acceleration timescales of the Kraichnan phenomenology is due to the fact 
that the overall wave intensity is given by 2it^, instead of 3u^ /2 as is for the Kolmogorov phenomenology. 

When Dkm ^ vf, corresponding to the fast diffusion limit, Tac — 'iD/v? and Tesc — {k^D)~^ . Then Tac/fesc — 



3D km/u ^ 1. The acceleration is negligible. With these asymptotic expressions for these timescales, the ratio of the 
acceleration to escape timescales decreases with the decrease of D and becomes close to unity near D ~ u/{3 



l/2l 



Since 
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D/c ( 1.0e15 cm ) D/c ( 1.0e15 cm ) 



Figure 3. The acceleration (solid) and escape (dashed) timescales as a function of the diffusion coefficient D for conditions at the 
transonic point xi in the downstream, where vp = u, for the Kraichnan (left with x = Hi ^ = 10.7 X lO^'^cm and Vy^ = 0.036(7) and 
Kolmogorov (right with x = 4.7, L = 8.0 X lO^^cm and Vj^ = O.OlSf/) phenomenology in Figure[T] U = 4000 km/s for both cases. The 
vertical line indicates the diffusion coefficient D = cl^/3. The thick lines are the exact results given by Equations II38I I and 1401 1. The 
thin solid line indicates the approximate acceleration timescale given by Equation l|f 6p . The thin dashed line indicates (k^D)~^ . For 
a magnetic field of _B = 14/xG, the diffusion coefficient for TeV electrons is about fO^^c cm. The c orresponding diffusion timescale over 
a length of ~ 3 X fO^® cm, the half- width of the variable X-ray filements observed with Chandra jUchivama et al~ll2007l ). is about one 
years. 



deviations from these asymptotic expressions occur near D ~ vp/km > u/km, efficient particle acceleration is only possible 
in the regime where D < vp/km- When Dkd <^ vp, corresponding to the slow diffusion limit, Tac — 3(5 — S)vp/[{S — 
S)Dkl-'kt'u^]. 



((5-3)u^ fDk, 



3D 



fDkrnV 
\ Vp J 



S~3 pk^D/vf 



dx- 



1 + 



3D 



dx- 



1 + x^ 



u-''D{Dkm/vpf^\{Al) 



where 3 < 5 < 5, which ensures the convergence of the integration. These results are in agreement with Figure [S] 

The acceleration and escape timescales are determined by D and the turbulence spectrum. Since we assume that D = Zdc/3 
for electrons with r-g < Id, the electron acceleration and escape timescales are independent of the energy for Vg < Id- When the 
ratio of acceleration to escape timescale is indepen dent of the par ticle energy, a power-law particle distribution is expected in 
the steady-state with the spectral index given by I Ptuskin 19881 ) 



P = 



/9 ^ Tac 
\ 4 Tesc 



(42) 



1/2 

Tesc I 2 

Since electrons with an energy less the proton rest energy can be scattered by whistler waves at small sc ales, the diffusion 
coefficient D for these particles can be much less than idc/3 and the corresponding spectral index p ~ 1 (jPetrosian fc Liu 
[2004 : Liu ct al. , 2006 ). The scatter will be dominated by the turbulent magnetic field before electrons reach the proton rest 
mass energy, the spectral break may appear below the proton rest mass energy. Without detailed treatment of these processes, 
we will assume that the electron distribution follows a power law with p = 1 for 7 < 10 in the following^ Above 7 = 10, 
the steady-state spectral index of the electron distribution is given by Equation 1)42^ since the diffusion coefficient D — ldc/3 
is assumed to be independent of the electron energy. For relativistic electrons with rg > Id, D > rgc/3. The increase of D 
with the electron energy E will lead to softening of the electron distribution toward higher energy dictated by the energy 
dependence of Tac/reac- The exact electron distribution can be obtained numerically. For the sake of simplicity, we assume 
that the electron distribution has a high-energy cutoff at Lorentz factor 



7c 



qBLv\ 



vaI min(ii_F, u) 
1 



for IK , 
for Kol . 



(43) 



where q is the elementary charge units. Then for a steady-s tate treatment, the di stribution of electrons escaping from the 
acceleration site may be approximated reasonably well with ijBecker. Le &: Dermer ii2006. : iPark fc Petrosion ,,1995 ) 



^ Assuming a spectral break at 7 = nip /me will lead to a spectral bump near this break energy, in conflict with radio observations. This 
suggests that, when the acceleration time is long and p is high, one needs better treatment of the electron acceleration at low energies. 
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Figure 4. Evolution of the acceleration efBcicncy rj (dotted), cutoff Lorentz factor 7c (dotted-dashed), spectral index p (dashed), and 
T = Tac/Tiifc (thin solid) in the downstream for the Kraichnan phenomenology with Vji = 0.036(7 and x = H (left) and Kolmogorov 
phenomenology with Vj^ = 0.018(7 and x = 4.7 (right). U = 4000 km/s. The particle acceleration is significant for r < 1. We only consider 
acceleration between the two vertical dashed lines indicating xi and X2. For 7c, we have assumed that B = 14^G, L = 10.7 X 10^^ cm 
(left) and B = 14/iG, L = 8.0 x lO^^^ cm (right). 
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Figure 5. Normalized steady-state nonthermal electron distribution f{x) produced at several locations in the downstream for the 
Kraichnan (left) and Kolmogorov (right) phenomenology in Figure |4] 



f{x,'y) cc 'y-^^^^e^p{-[Tac{D)/r.sc{D)f^}/r.sc{D) , 



(44) 



where 



D{'y,x) = 



dd{x)/3 
crg(7)/3 : 



: 7mccV(3gB) 



for 7 < 7c (a;) , 
for 7 > 7c (a::) ■ 



(45) 



We note that the shape of the distributio n function near 7c can also be affected by the energy loss processes (jStawarz fc Petrosian 



2008 



1 IVannoni et al. II2OO9I: iBlasi Iboioh. In our model, it is mostly det e rmined by the balance between the acceleration and 



escape processes (jPark fc Petrosion Ill99i : iBecker. Le fc Dermer I i2006l : IZirakashvili fc Aharonian 1120071 ). The shape of this 
high-energy cutoff has significantly effect on the fitting parameter, especially the value of 7c. We are carrying out detailed 
numerical investigation of the particle distribution and the results will be reported in a separate publication. 

Although the parallel propagating waves in the dissipation range may not contribute to the particle scatter significantly, 
energies carried by these waves are only accessible to relativistic electrons. At very small scales, these waves resonate with 
thermal background electrons giving rise to a preferential acceleration of electrons. One therefore may assume that the ratio 
of the dissipated energy carried by non-thermal electrons to that of the thermal ions is proportional to the ratio of the energy 
density of parallel propagating fast-mode waves to that of the magnetic field: 
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Figure 6. Distributions of nonthermal electrons Fix) at several locations in the downstream for the Kraichnan (left) and Kolmogorov 
(right) phenomenology in Figure |4] 



_ eg (fed) _ f e''''v\/(2^f/\sVF for IK, 

2 "\ e5/%^/(27r)^/2^;5 for Kol . ^ ' 

where e ~ 2.72 is the base of the n atural logarithm . A more quantitatively treatment of this issue may address the electron 



injection processes self-consistently (jEichler Ill979bl ) 



To have efficient acceleration of relativistic electrons, the turbulence decay time 

„r 1^ f tci&-k{u,vf) lu for IK , 
1 ' forKoi. 

and the remnant lifetime Tnfe should be longer than the acceleration time. As we will see below, the turbulence decay time 
is always longer than Tufe in the subsonic phase for SNR RX J1713. 7-3946. There are at most two locations xo and X2 with 
xo < X2 in the downstream, where r = Tac/Tiifo = 1- Figure |4] shows the evolution of rj, 'jc, p, Tac/Td and r = rac/Tiifo in the 
downstream for U = 4000 km/s. The profiles of vf /U and u/U only depend on va/U. So is the profile of ry. The profiles of 
T and p also depend on the absolute value of U. To obtain 7c, one needs to know L and B as well. Most SA occurs near the 
sonic point xi, where vf = u. In the late subsonic phase, u <^ vf, the SA is insignificant since most of the free energy of the 
system has been converted into heat. The characteristic length of the magnetic field is also long far downstream due to the 
weak turbulence, which implies long electron scatter and acceleration timescales. 
Then the distribution of non-thermal electrons in the downstream 

F{x,'y)= [ /(x',7)r?(x')(4Q/mec'(7)da;' (48) 



where 7/(k', 7)d7 — 1, and jmec'^F{x,j)d'y gives the energy density of non-thermal particles at x. Figures [5] and 
[6] show the normalized electron distribution / and F for parameters in Figure |4] at several locations in the downstream, 
respectively. We note that the electron distribution has a rather gradual high-energy cutoff due to the weak dependence of the 
acceleration and escape timescales on the spatial diffusion coefficient D, which is proportional to the electron energy above 
the cutoff energy 7cmeC^. This gradual cutoff results in a broad emission component due to inverse Comptonization of the 
low-energy background photons by high-energy electrons, which can fit the observed broad TeV emission spectrum from a 
few SNRs. 



5 APPLICATION TO SNR RX J1713.7-3946 AND TIME-DEPENDENT MODELS 

Here we use the SNR RX J1713. 7-3946 as an example to demonstrate how the SA by fast-mode waves accounts for the 
observed broadband spectrum. The electron distribution produced by the SA is given by equation (|48|) . where the integration 
over x' should stop at X2- Given the evolution history of the SNR, the nonthermal electron distribution will vary along the 



radia l direction. A detailed modeling o f the explosion is necessary to take this effect into account properly (jCowsik fc Sarkar 



198 



3; Zirakashvili fc Aharonian 2O10l l. Here we treat the volume of the emission region Ve as a free parameter to control 



the normalization of the emission spectrum, which is appropriate for SNRs where the nonthermal particles appear to be 
concentrated near the SF. By adjusting U, B, va, X) L, and 14, one can use the corresponding electron distribution F(x2) 
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Figure 7. Best fits to the observed spectrum of tlie SNR RX J1713. 7-3946. Tlie X-ray data points are obtained from iTa naka et al. I 
l|2008h . The other data points and the sensitivity limits of different high-energy telescopes are the same as those in lLiu et al. I l l2008al '). 
The left and right panels are for the Kraichnan and Kolmogorov phenomenology in Figure [J] respectively. The dashed line is for a 
simple powrer-law model with a gradual high-energy cutoff F('y) oc ^^-^^ exp(7/7c)"'^/^, where 'fcirieC^ = 3.68 TeV. The solid lines are for 
the fiducial models. The low and high e nergy spectral peaks are produced through the synchrotron and inverse Compon scatter of the 
background photons l lPorter et al. |[2006t) . respectively. The preliminary data fr om Fermi w ere not considered in the fit and are included 
here to show the agreement between model predictions and this observation l lFunkll2009l) . The Fermi observation seems to favor the 
Kraichnan model. 



Table 1. Model Parameters 



model 


U 


B 


va/U 


X 


L 


v,/v 




P 




(km s^-"-) 


(mG) 






lOi'^(cm) 




10'''^{erg) 


10-2^{g cm-3) 


IK steady 


4000 


14.0 


0.036 


11 


10.7 


0.42 


7.69 


7.52 


KOL steady 


4000 


14.0 


0.018 


4.7 


8.0 


0.016 


6.96 


30.1 



to fit the observed spectrum of the SNR RX J1713. 7-3946. Figure [7] shows the best fit with the corresponding parameters 
listed in Table [TJ where V is the enclosed volume of the SNR SF. Comparing to the thin-dashed line, which is derived by 
assuming an electron distribution oc exp — 

(7/7c)''", the high -energy cutoff of the Kraichnan model is more gradual than 
the Kolmogorov model, which makes the former spectrum broader. The emission volume is 0.42 (0.016) times the volume 
of the SNR for the Kraichnan (Kolmogorov) phenomenology, which is compatible with observations. The values of x ^^re 
greater than 1, which suggests that the particle spatial diffusion coefficient might be enhanced significantly by incompressible 
turbulence motions. The total energies of nonthermal electrons are also comparable to the magnetic field energy for a uniform 
magnetic field within the remnant, suggesting near energy equipartition between the nonthermal electrons and the magnetic 
field. 

There are six parameters in the model: B, U, va, X) L, and 14. The observed radio to X-ray spectral index. X-ray to 
TeV fiux ratio, location of the X-ray cutoff, and bolometric luminosity of the source give four constraints, which leads to 
two more degree of freedom. Our model fit to the spectrum therefore is not unique, x is determined by the coupling between 
incompressible and compressible turbulence motions. B is well constrained by the ratio of the X-ray to TeV fiux. To reproduce 
the observed spectral shape, the profiles of p, 7c, and rj should not change significantly, which implies that v\c? oc u^" and 
L oc u*' /v\ at the transonic point for the Kraichnan phenomenology. For the Kolmogorov model, v\c? oc and L oc /v\. 
For Va ^ U , u is proportional to U, we find that nearly identical emission spectra can be obtained by adjusting U and va 



I Liu et al. Il2008bl ). Since the turbulence decay time is longer tha n the supernova lifetime, the steady-state treatment is also 



justified. A time-dependent treatment gives identical parameters ( Becker. Le fc Dermer 20061 ) 



6 DENSITY AND TURBULENCE GENERATION 



The primary discrepancy betwee n these models and th e obse rvations are the relatively high densities of the downstream 
plasma. From X-ray observations, Cassam-Chenai' et al. ( 20041 ) inferred an upper limit for the electron density of 0.02 cm~^. 
The corresponding mass density is about 3.3 x 10~^® g cm~^, which is comparable to the densities inferred with the Kraichnan 
phenomenology but lower than those inferred with the Kolmogorov phenomenology. On the other hand, the electron tempera- 
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ture c ould be much lower than the ion temperature in the shock downstream ( Zirakashvili fc AharoniajTlboiol ). Morlino et al 



1 2OO9I1 and lFang etHI (|2009l ) lave argued that the density in the downstream can be as high as 0.5 cm , corresponding to 



a mass density of 8.4 x 10 ^^g cm ^. 

Given the age of Tufc = 1600 years and a radius of 7? = 10 pc at a dist ance oi D ^ 1 kp c, the corresponding average 



speed of the shock front is 6100 km s ^. With the self-similar solution of I Chevalier 

I l|l982l ). we infer a shock speed of 



[(n — 3)/(n — s)]6100 km s ^, where n > 5 and < s < 3 are the power-law exponents of the density pro file for the ejecta 



and the ambient medium, respectively. Observations give an upper limit of 4500 km s ^for the shock speed ( Uchivama et al 



20071 ). From the self-similar solution, the shock speed must be higher than 2400 km s ^ (for s — Q and n = 5). If ions are 



preferentially heated by the shock, the corresponding ion temperature Ti will be higher than Zm-pU^ /l&kB, > 1. 3 x 10^ K. The 



electron temperature should be higher than that given through Coulomb coUisional energy exchange with ions (jHughes et al. 
2000h : 



Te > 2.1 X 10''(Tiife/1600yr)^''^(ne/cm-^)^/^(r,/1.3 x 10**K)^/'^K, (49) 

where Ue is the electron number density. The corresponding bremsstrahlung luminosity is £ > 5.2 x 10'^''(ne/0.5cm~^)^^''^ erg 
s~^, which is comparable to the l uminosity of the obse rved nonthermal X-ray emission|f] We therefore expect strong thermal 
emission with such a high density. Morlino et al. ( 20091 ) obtained a very low thermal bremsstrahlung luminosity by arbitrarily 



adopting an electron temperature 100 times lower than the ion temperature. As shown above, considering the electron ion 
Coulomb coUisional energy exchange, the electron temperature will not be that low and significant thermal X-ray is expected 



with a density of 0.5 cm , except that cooling of the shock front by cosmic ray ions dominates (|Zirakashvili fc Aharonian 



2OI0I) . The highest electron density given by our models is about 0.2 cm . The corresponding thermal X-ray luminosity 



will be reduced by nearly one order of magnitude and should be in agreement with observations. Detailed modeling of the 
supernova explosion and the thermal emission is needed to see the validity of these models. 

The model inferred density may also be reduced by considering the acceleration of electrons by large-scale structures 
in the downstream and the acceleration in the supersonic phase, where the first-order Fermi acceleration is also possible. In 
this paper, we consider the electron acceleration by the fully developed turbulence in the subsonic phase. It assumes that 
once the electrons diffuse over a scale of the turbulence generation length L, the acceleration stops. As shown above, the 
turbulence evolves in the downstream. In a more self-consistent treatment, one may use the turbulence properties to derive 
nonthermal electrons injected into the downstream flow by small scale plasma waves and consider the further acceleration of 
these electrons as they diffuse spatially in the downstream. The scatter mean-free-path of these particles are determined by 
the properties of turbulence. The electron acceleration stops only after they diffuse into upstream or far downstream, where 
the turbulence becomes insignificant. If these effects lead to a harder overall electron distribution, the Alfven speed needs to 
be increased to fit the observations, leading to a lower density. 

From these models studied here, we see that, to have efficient SA, both high-speed waves {vf — u) and short scatter 
mean-free-path are required. Quantitatively, one needs c?l^ / L^v? to be on the order of unity so that the acceleration and 
escape time scales of relativistic particles are comparable, u is constrained by the shock speed. A short scatter mean-free-path 
is achieved by the reduction of the characteristic length of the magnetic field, which also determines the maximum energy of 
the accelerated particles. In these models, turbulence motions are invoked to reduce the characteristic length of the magnetic 
field. It is obvious, such a mechanism is only possible for strong turbulence, where the turbulence speed is higher than the 
Alfven speed. The turbulence speed is determined by the shock speed u < (3/16)^^^(7, which is less than 1949 km s~^ for 
SNR RX J1713. 7-394. Therefore va < 1949 km s"^ and we obtain a low limit for the mass density from the inferred magnetic 
field of 14 fiG: p = B^/Attva > 4.1 x lO'^^g cm"^(B/14^G)^((7/4500km s"^)"^, which corresponds to an electron density of 
~ 0.0002 cm-^(B/14^G)^((7/4500km s~^)"^ 

Therefore further reduction of the density can be achieved by considering the generation of the turbulence and its effect 
on the turbulence spectrum, i.e., the dependence of the eddy velocity on the spatial scale. Indeed, in the models considered 
above, we assume that the turbulence is generated in a very narrow spatial range instantaneously at the SF and an inertial 
range develops. Since the large-scale eddy speed is comparable to the bulk velocity of the downstream flow moving away from 
the shock front, the region with x < 0.5 should be considered as the turbulence generation phase. This is an intrinsic limitation 
of the above treatments, which focus on the averaged properties of the downstream flow without addressing the turbulence 
generation process. Significant particle acceleration occurs within x < 0.5, i.e., the turbulence generation phase. It is also 
possible that the turbulence is generated over a broad spatial range and/or the turbulence is not isotropic at large scales. One 
then expects a turbulence spectrum shallower than the initial range spectrum. For example, for a turbulence spectrum of 

W = [3u/8Tv{g - 3)](27r/L)»-^fc-«-'' (50) 



Here we have assumed an emission volume one quarter of the volume enclosed by the remnant shock front. 



Stochastic Electron Acceleration in TeV SNRs 15 



with 1 < g < 1.5, □ where the normahzation is chosen so that the turbulence energy density is given by {3/2)pu^, the eddy 
speed can be redefined as 

= {[STvig - l)/3]k'W}'^^ = ^(fcL/2^)(^-«)/2 . (51) 

The eddy speed is comparable to the Alfven speed at the characteristic length of the magnet field I, we then have 12c^/^/ujL^ — 
\2<?v'^l'-''-^'^ lv)u^l^^~^^ ~ 1. Therefore vaIu ~ (v; jVlcf^'^^l'^ . Since ~ u ~ [/, we have vaIU ~ ([//c)'^-^'/^ The Alfven 
speed can be comparable to the turbulence speed for a shallow turbulence spectrum with g approaching 1. Thus stochastic 
electron acceleration can account for observations of SNR RX J1713. 7-394 as far as the mass density of the shocked plasma 
is greater t han 4.1 x 10~^^g cm~'^. Detailed studies of the turbulence generation and the associated pa rticle acceleration are 



warranted (|Lucek fc Bell Iboool : lllededal et ^1120041 : [ciacalone fc Jokiph II2OO7I : iNishikawa et al. II2OO9I ') 



7 SUMMARY AND CONCLUSIONS 

In the paper, we study the SA of electrons by a decaying turbulence as produced and/or enhanced by strong non-relativistic 
shocks and carried away from the SF with the downstream flow. It is shown that, to have signiflcant particle acceleration, the 
turbulence must cover a large spatial scale so that the particle acceleration time may be shorter than the turbulence decay 
time. To account for observations of a few STTSNRs with the leptonic scenario for the TeV emission, fast-mode waves need 
to be excited in the subsonic phase. Given the turbulent nature of the downstream flow, fast-mode waves may prevail in the 
downstream. We show that the SA by large-scale acoustic (fast-mode) waves can account for the observations. 

There are four basic model parameters, namely, the magnetic field, mass density, shock speed, and the turbulence gener- 
ation scale. O bservations of a few SNRs and radio galaxies have shown that t he particle acceleration may change dramatically 
along the SF jRothenflug et al. II2OO4I : iRevnolds II2OO9I : ICroston et al. II2OO9I ). This variation has been attributed to a large- 



scale magnetic field in the DSA model as the acceleration efficiency varies with the angle between the magnetic field and the 
shock normal. With our model, this variation is likely caused by a quite different mechanism. Detailed comparative studies 
should be able to distinguish these models. 

The particle acceleration is very sensitive to the magnetic field. To produce nonthermal particle distribution in compatible 
with observations, c^w^/it^" and c^v^/u^, where u is the large-scale eddy speed near the shock front, should be on the order 
of 1 for the Kraichnan and Kolmogorov phenomenology, respectively. The high-energy cutoff of the particle distribution is 
determined by the magnetic field and the turbulence generation scale. Weaker fields will lead to lower cutoff energies and 
harder spectra. The particle acceleration may be turned off completely for strong fields due to the increase of the characteristic 
length of the magnetic field, and therefore the particle scatter mean-free-path. If the magnetic field is predominantly generated 
by the streamline of nonthermal particles upstream, the models then imply va ~ U^^'^ and va ~ V^^^ for this dynamo process 
and for the Kraichnan and Kolmogorov phenomenology in the downstream turbulence, respectively. 

Assuming that the turbulence is isotropic and generated in a narrow spatial scale, the model inferred densities of the 
downstream fiow may be so high that thermal X-ray emission becomes observable, in confiict with observations. Although 
the thermal X-ray emission may be suppressed by a lower electron temperature due to dominance of the cooling by cosmic 
ray ions, we find that a low density is also possible if the turbulence is not isotropic or generated over a broad spatial scale so 
that the eddy speed has very weak dependence on the sp atial scales. Detailed mo deling of the progenitor and the evolution 
history of the remnant may help to constrain the density ( Cowsik fc Sarkar 1984h . 



The application of the models to the SNR RX J1713. 7-3946 also suggests energy equipartition between the magnetic 
field and the acceleration electrons. If the rest of the shock energy is dissipated as heat in the downstream, then the overall 
electron acceleration efficiency will be ~ va/u^, which is inversely proportional to the plasma /3 of the downstream flow. For 
the Kraichnan phenomenology, v^/u^ ~ {u/c)^'''^ and for the Kolmogorov phenomenology v^/u^ ~ (m/c)^^^, acceleration is 
more efficient for stronger shocks, which also produce hotter downstream plasma. However, the dependence of the acceleration 
efficiency on the shock speed is rather weak. These ma y have significant implications on the origin of cosmic rays and their 
connection to the properties of the interstellar medium I Reynolds 20081 ). 



With the fast-mode wave turbulence model studied in this paper, relativistic electrons may also be accelerated through 
the first-order Fermi mechanism as in the DSA model, especially in the supersonic phase. The high-energy cutoff can still result 
from decoupling of higher energy particles with the background magnetic field as their gyro-radius exceeds the characteristic 
length of the magnetic field. This is quite different from the DSA mode ls, where the high-energy cutoff is related to a finite 



lifetime of the shock, shock curvature, or efflcient energy loss processes I Zirakashvili fc Aharonian 201Cll ) 



^ For g < 1, the eddy speed increases with the decrease of the spatial scale and the turbulence energy is not dominated by the large 
scale eddies. One then needs to introduce a spectral break at certain small scale. Such a complex scenario is not well justified from both 
observational and theoretical point of views. 
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